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ABSTRACT 

We study here the random diffusion model. This is a continuum model for a conserved 
scalar density field (p driven by diffusive dynamics. The interesting feature of the dynamics is 
that the hare diffusion coefficient D is density dependent. In the simplest case D = D + Di6(j) 
where D is the constant average diffusion constant. In the case where the driving effective 
Hamiltonian is quadratic the model can be treated using perturbation theory in terms of 
the single nonlinear coupling Di. We develop perturbation theory to fourth order in Di. 

The are two ways of analyzing this perturbation theory. In one approach, developed by 
Kawasaki, at one-loop order one finds mode coupling theory with an ergodic-nonergodic 
transition. An alternative more direct interpretation at one-loop order leads to a slowing 
down as the nonlinear coupling increases. Eventually one hits a critical coupling where 
the time decay becomes algebraic. Near this critical coupling a weak peak develops at a 
wavenumber well above the peak at q — associated with the conservation law. The width 
of this peak in Fourier space decreases with time and can be identified with a characteristic 
kinetic length which grows with a power law in time. For stronger coupling the system 
becomes metastable and then unstable. 

At two-loop order it is shown that the ergodic-nonergodic transition is not supported. It 
is demonstrated that the critical properties of the direct approach survive going to higher 
order in perturbation theory. 
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I. INTRODUCTION 

We study here a dynamical system, the random diffusion model (RDM), under going 
diffusive dynamics with a field dependent diffusion coefficient. This model serves as a very 
simple model for the dynamics of the density field in colloidal systems. The motivation 
for studying this model comes from facilitated spin models where the kinetic coefficient 
is density dependent and leads to significant slowing down for dense systems. Here the 
expectation is that a properly chosen bare diffusion coefficient leads to significant slowing 
down as the density increases. 

There has been much speculation but relatively few solid results in establishing the ex- 
istence of a mode coupling theory (MCT)P ergodic-nonergodic (ENE ) transition in field 
theoretic models of the liquid-glass transition. The RDM is a candidate for the simplest 
such model. At one-loop ordei^ the self-consistent theory can be organized in two ways. 
One approach, involving a rearrangement due to KawasakP, leads to a conventional ENE 
transition for a large enough nonlinear coupling. A second, more direct approach, leads, on 
increasing the same dimensionless coupling, to slowing down of the system for wavenumbers 
well away from zero. One eventually reaches a coupling where the system produces a peaked 
dynamic structure factor at wavenumbers well away from zero. We can call this a prepeak 
since we expect it to show up at wavenumbers below those characterizing the first peak in 
the static structure factor. The width of this peak corresponds to a kinetic length which 
increases algebraically with time. 

The key point here is that the RDM is simple enough that the associated perturbation 
theory, in terms of a single dimensionless expansion parameter, can be worked out at two- 
loop order. We find that the ENE transition does not survive in the two-loop theory for 
reasons that may be generic. In a direct self-consistent treatment of the model we find that 
the new peaked state is maintained at two-loop order. 

We focus on field theoretical models for the dynamics of dense fluids since they offer 
the best hope of a self-consistent theory. This hope includes the possibilities of higher-order 
computation and the determination of four-point correlation functions. As mentioned above, 
the computation at higher order is necessary to establish the stability of any ENE transition 
found at one-loop order. We also want to compute multi-point correlation functions since 
there is speculatiorP that they offer information about a growing length as one approaches 
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the ENE transition. 

Despite a number of papers (see below) discussing mode coupling theory (MCT) from 
the point of view of field theoretical models, the situation is unclear. We do not really 
know which models have a transition and which do not. There is worl<P suggesting that 
nonlinear fluctuating hydrodynamics offers viable kinetic models for studying the dynamics 
of dense fluids and can lead to the ENE transition. Das and Mazenko (DM)^ introduced a 
field theoretical model with density and momentum fields. They showed from general non- 
perturbative considerations and a one loop calculation that the conventional mode coupling 
transition is cutoff. Schmitz et aP found a cutoff in a slightly simpler model. Gates and 
Ramaswamy^, using heuristic reasoning, argue that these cutoffs are not effective in the DM 
model. 

A set of slightly simpler models (involving only the density field) were introduced by 
DearP and KawasakP^ to describe the overdamped diffusive dynamics in colloidal systems. 
Miyasaki and Reichmari^ studied the Dean-Kawasaki (DK) model using the MSR methocP. 
They found a nonlinear fluctuation dissipation theorem (FDT) connecting propagators and 
correlation functions which made even the one loop theory difficult to interpret. Things are 
complicated by the use of the MSR method which requires field doubling in carrying out the 
perturbation theory. Andreanov, Biroli, and Lefevre (ABL)— document that nonlinear terms 
in the effective hamiltonian generate a nonlinear FDT and make systematic perturbation 
theory very difficult. They suggest introducing auxiliary fields to solve this problem but 
were unable to construct a sensible one loop approximation. Kawasaki and Kim-^, taking 
a similar approach, were able to find a one-loop approximation which does lead to a ENE 
transition. 

Given the uncertainty in the analysis of the DM and DK models (Do they have an ENE 
transition?), we need to analyze a simpler model which does have an ENE transition at one 
loop order. The idea is to check whether this solution is stable at two-loop order. These 
technical problems suggest that one needs to study a model that has a linear FDT. 

The RDM is related to the DK model. It is the simplest nontrivial realization of the 
hindered diffusion model^ introduced earlier. The physical motivation for this model is from 
facilitated spin model ^ * ^^ * ^^ * ^^ * ^*^ where the kinetic coefficient in a lattice model dynamics 
depends on the local environment in a constraining manner. In a continuum model, with 
a conserved density, the analog is a density dependent diffusion coefficient. In both models 
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one can have strong kinetic slowing down despite having trivial, "noninteracting" static 
equilibrium behavior. 

The random diffusion model model has a single identifiable small parameter. As discussed 
in referencdlSl the source of nonlinearities are in the density dependence of the bare diffusion 
coefficient. In the simplest case the bare diffusion coefficient is of the form 

D{<P)=Do + D^(t> (1) 

and the perturbation theory is in powers of Di. In the simplifying case where we assume 
the static structure, in our coarse grained system, is a constant up to a cutoff A, called here 
the structureless approximation, we find that the dimensionless coupling constant is given 
by 

where 0o = (0) is the average density and 

S = {{5ct>f) (3) 

is the local fluctuation in the densitj^. 

The RDM shows an ENE transition at one loop {g^) order but the transition appears 
to be inconsistent with the theory at two-loop order. There is an alternate more direct 
approach to the perturbative treatment of the RDM. In this model, as a function of increasing 
coupling, (yf, one finds a slowing down. For coupling g < g*, where g* is the critical coupling, 
there is a cross over from exponential to algebraic time decay for a band of wave numbers 
away from zero wavenumber. Indeed certain wavenumber components decay to zero more 
slowly than others and a small amplitude peak develops in the dynamic structure factor. 
This structural peak has the form 



The width of this small amplitude peak decreases with time thus giving a length v B which 
increases algebraically with time. The amplitude A decreases with time and, after a brief 
initial transient, Qo, is time-independent. A shows power-law behavior in time for g near 

g*- 

For g > g* the system is slow but eventually unstable. The small peak contribution, for 
long enough time, begins to grow and the system eventually blows up. It is not unreasonable 
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to assume that the unstable system represents the nucleating solid phase. The model must 
be extended with the appropriate static behavior if one is to stabilize the nucleated solid 
phase. 

It has been traditional to use the MSR method to develop perturbation theory for dy- 
namical models such as the RDM studied here. This method has the distinct advantage 
that perturbation theory can be developed in terms of the physical correlation and response 
functions. In the RDM the correlation and response functions are linearly related and the 
calculation at one-loop order is manageable. The calculation at two-loop order is extremely 
complicated by sums over the labels differentiating fields from response fields. The Fokker- 
Planck description has the advantage that the bare perturbation theory expansion is formally 
transparent, the static behavior is easy to sort out, and one does not have the frequency 
integrals found in the MSR method. The disadvantage is that one has to replace the bare 
correlation functions by their full counterparts by hand. One is helped by the knowledge 
from the MSR approach that such a renormalization (resummation) exists. 



II. RANDOM DIFFUSION MODEL 



We discuss our model in the context of a Fokker-Planck description. The equihbrium 
intermediate dynamic structure factor is given by 



C(qi,q2;t) = / P(</.)iy^</)(q2)e-^^V(qi) 



(0(q2)e-^^V(qi)) 



(4) 



where 0(qi) is the Fourier transform of the fundamental field 5(j) in the theory, the equilib 
rium probability distribution is given by 

z 



(5) 



where the effective Hamiltonian Ti^ can be taken to be quadratic in 0: 

'H^ = j c/''xirf''x2^50(Xi)x-l(Xi -X2)50(X2) (6) 

and 50(xi) = 0(xi) — 0o- The adjoint Fokker-Planck operator for our model is given by 



50(x) 50(x) 



r0(x,y) 



50(y) 



(7) 



wherd^Sl 



r^(x,y) = V,-V, (D(0(x))5(x-y)) 



(8) 



and the bare- diffusion coefficient is taken to be of the simplest nontrivial form 



L'(0(x)) = Do + A</'(x) 



(9) 



A more complicated and physical form for -D(0) was studied in RefJTSl 

Our model can also be written as a field theory of the MSR!^ type. The MSR action is 
given in this case by 



A = J d'^x dt 



- V. D(0)Vi 



(10) 



where is the MSR auxiliary response field. 
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III. MEMORY FUNCTION FORMALISM 

We use here a memory function formalism in the Fokker-Planck description. This ap- 
proach was first fully developed in RefJ251 for kinetic theory and later appliecP^ to the fluc- 
tuating nonlinear hydrodynamics of smectic A liquid crystals. A significant advantage of 
the method is that it allows one to treat interactions expressed in terms of static averages. 
Thus, in the present problem, it is static averages not the bare diffusion coefficient which 
appears in the theory. The structure of this type of theory was investigated in some detail 
by AndersorpSl. 

Let us work with the Fourier-Laplace transformed time correlation function 

oo 



C(qi, qi; 2;) = -i / dte''^C{qi, qi; t) 
Jo 



= (0(q2)i?(^)0(qi)) , (11) 
where the resolvant operator is given by 

/■oo - 

R{z) = -i / dte'^'^'^^^' = [z + iDJ-^ . (12) 
Jo 

Using the identity 

zR{z) = 1 - R{z)iD^ (13) 



in Eq.(ll), leads to the kinetic equation 



^^^K(qi,ki;2;)C(ki,q2;z) = C(qi,q2) . (14) 



The memory function, K, is given by 

r(qi, q^; z) = (27r)'^5(qi - q2)r(qi; z) = J (^^(Qi, ki; z)C{k^, 



r(^)(qi,q2) + rW(qi,q2;z) . (15) 



The static part of the memory function is given by 

r(^Hqi,q2) = (0(q2)*^00(qi)) = 



(0(q2)t;(qi)) 



(16) 



where the current v is defined by 

v{qi) = iD^(f){qi) , (17) 
and the dynamic part of the memory function is given by 

r^'^^(qi,q2;2;) = r(qi,q2;2;) + r^„5(qi, qa; z) (18) 

where 

T{q,,q2;z) = -{v{q2)R{z)v{q,)) (19) 
and the subtraction part is given by 

r.„.(qi,q2;^) = /^/ ^,w{q2M,^)c-\k2M;^)wiqM (20) 

where 

W^(qi,ki;^) = (0(ki)/?(^)t;(qi)) . (21) 
Using standard arguments we can show that the physical diffusion coefficient is given by: 

D„ = hm hm -i^T(q, z) . (22) 

IV. BARE PERTURBATION THEORY 

A. Two-time quantities 

In this section we show how to set up perturbation theory for the dynamic structure 
factor 

C(qi,q2;^) = (0(q2)/?(^)0(qi)) (23) 



where the resolvant operator is defined by Eq.(12). We are then interested in carrying out 



perturbation theory where the FP operator can be written as the sum 

D, = Df + Df , (24) 
where the zeroth-order contribution is given by Eq.(??) with replaced by 

Tf{x,y) = DV.-Vy5{^-y) (25) 



and the interacting contribution is given by Eq.([7]) with F,^ replaced by AF^ given by 

AF^(x, y) = V, ■ V, (Di50(x)(5 (x - y)) . (26) 
We then use the operator identity 

R{z) = [z + zD'^^ + zDf]~' 

= [z + zDi°)]- + [. + ^DfY\-^Di'')[z + ^Df + zD^'Y' 

= Mz) + Ro{z){-tDf)R{z) (27) 
which defines the zeroth-order resovent 

Roiz) = [z + tD^;^Y' ■ (28) 
Using this result in the correlation function and iterating gives 

R{z) = Ro{z) - R,{z)iDfRo{z) + R,{z)iD^^ R^{z)iD^'^ Ro{z) + ... 



oo 



= i?o(^)EH^r^o(^)]" . (29) 

n=0 

For correlation functions we have the expansion 

C^(qi,q2;^) = (0(q2)i?(^)0(qi)) 



oo 



EC('^)(qi,q2;^) (30) 



n=0 

where order by order 



C(°)(qi,q2;^) = (0(q2))i?o(^)0(qi)) (31) 

C(i)(qi,q2;^) = (0(q2)i?o(^) H^r^o(^)]0(qi)) (32) 

C(2)(qi,q2;;^) = (0(q2)i?o(^)H^i'^i?o(^)]'</'(qi)) (33) 

C(")(qi,q2;z) = (</>(q2)i?o(^)[-^^^'^i^o(^)]>(ql)) • (34) 



The first step in the analysis is evaluate Ro{z)(j){cii). In Appendix A we show 

Ro{z)(f){qi) = To{qi, z)(f){qi) 

where 

To{qi,z) = [z + iLoiqi)]'^ 



and Lo{q) is given by Eq.(84) below. 



On the left in each expression for C*-"^ we have 

{4>{q2)Ro{z) ...) = {[Roiz)<Piq2)] • • •) = T.iq^, z)(0(g2) • • •) • 
At the various orders we have 

C^^°n?i,92),^)=To(gi,z)(0fe)0(gi)) 

C^'\qi,q2),z)=Toiq2,z){<P{q2)[-tDi;^]<Piqi))To{quz) 

C^^\q,,q,),z)=Uq2,zmq2)[-iD^l^Ro{z)t~^^^^ 
The zeroth order solution is explicit after identifying 

C^'^{qi,q2)) = (0(g2)0(gi)) = A;BTx(gi)(2vr)^(5(qi + q2) . 
If we introduce the interaction part of the current, 

v^'\q,) = tD^'Uiqi) , 
then for the higher-order contributions 

C^'\q,,q,),z) = -To{q2,z){<Piq2)v^'\q,))To{q,,z) 
C^'\qi,q2),z) = Uq2,z){v^'\q2)Roiz)v^'\qi))Uq„z) 

and 



Let us look at the nonlinear term v^^\qi). In coordinate space: 

6 



J 50(X2) 



where AF^ is given by Eq.(26) and Ti^ by Eq.([6]). Inserting these expressions into 
leads to the cubic vertex 



t;(^)(xO = -^Ev:,(A 



lO</)(Xil</)alXi 



where 



1^{^,) = V:J rf'^a;3X"'(x2-X3)(50(x3) . 



Taking the Fourier transform gives the cubic interaction: 

-''\^) = f^f ^F(q,k2,k3)0(k2)0(k3) 



where 



{2ttY J (27r) 



V{q, k2, h) = ^Diq ■ A{k2, h)i27iy6{q - k2 - k3) 



It will also be useful to write the vertex in the alternative form 



V{q, k2, ks) = V{q, fc2)(27r)'^5(q - k2 - k3 



and 



where 



V{q, k2) = 2^iq • {k2X-\k2) + (q - k2)x-'(q - ^2] 



Since v^^\q) is even in 50 we have 



and from Eq. ( 43 ) 



C^'\qi,q2,z) = 0. 
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Turning to Eq.(44) we have at second order in the couphng: 

d'^h d'^k2 d'^h d'^ki 



C^'\q,,q2,z)=To{q2,z) J 



{27rY {271^ {2ttY {27rY 



X V{q2, ks, k^)M^^\ki, k2, /C3, ^4; z)V{qi, k,, k2)To{qi, . 



(56) 



and 



M^'\k,,k2,ks,k,-Z) = (0(k3)0(k4)i?o(^)0(ki)0(k2)) 



(57) 



At higher order we have generally 



C^^\q,,q2,z)=To{q2,z) 



d'^ki d'^k2 d'^ks d'^k^ 
(2^(2^(2^(2^ 



X V{q2. ks, ki)M^'''\ki, k2, ks, ki- z)V{qi, k2)To{qi, z) 



(58) 



where 



M^''\kuk2,ks,k,;z) = (0(k3)(^(k4)i?o(^)H^f i?o(^)]^"-'V(ki)0(k2)) . (59) 



To go further with M^"^^ we need the result from Appendix A 



i?o(^)0(ki)0(k2) = To(ki, ka; z) [0(ki)0(k2) - ^(ki, ks)] 



C'(ki,k2) 



(60) 



Using this result in M*^^)(ki, k2, ka, k4; 2), Eq.(57), we can then do the static average over 
gaussian fields to obtain 

M(2)(ki,k2,k3,k4;;2) = ro(ki,k2;z) 



[^(ki, k3)C'(k2, k4) + ^(ki, k4)C'(k2, ka)] 



+ C'(k3,k4) 



C'(ki,k2 



Putting this result back into Eq.(56) and using the result 



(61) 



y(q, ki, k2)C'(ki, ks) ^ qJiq - ki - k2)S{k, + k2) = 



(62) 



rives 

C(^)(q„ q„ .) = To(q„ .) / ^^^^^(^^^ ks, 1.4) 



X ro(ki, ks; z)2C'(ki, k3)(7(k2, k,)Viq,, h, h)n{q,, z) . 

Using the 5-functions in the vertices and static correlation functions allows one to 
of the k integrals and obtain 

C(2)(gi,g2,z) = 2(27r)'^5(qi + q2)To(-gi,^) j ^l^(-qi, -ki)To(ki, qi - ki 

xC'(ki)(7(qi-ki)V(qi,ki)ro(gi,^) . 

We return to this expression below. 

Going to the third-order contribution we must evaluate 

M^''\k^,k2,h,k^;z) = (0(k3)0(k4)i?o(^)H^i'^i?o(^)]0(kl)0(k2)) . 

since ID'^^^ is odd in 50 it is easy to see that 

M^^\ki,k2,k3,ki;z) = 

and 

C^'\ki,k2,ks,k^;z)=(} . 
At fourth order we must evaluate 

M^'\k,,k2,k3,k,;z) = (0(k3)0(k4)i?oWH^i'^i?oW]V(ki)0(k2)) 

= ([2D('^i?o(^)0(k3)0(k4)]i?o(^)[^^r^o(^)</'(kl)</.(k2)]) . 



We find immediately, using Eq.(60), that 

zDl'^Ro{z)4>{k,)4>{k2) = To(ki, k2; ^)zDj,'V(ki)0(k2) . 
We have from appendix B 

^Dj,'V(ki)0(k2) = ^7(^)(kO0(k2) + t;(^)(k2)0(ki) + 5(ki,k2) 
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where 



^(ki, ka) = 2ip~^Diki ■ k20(ki + k2 



Using Eq.(70) twice in Eq.(68) gives 



M(^)(fci, h, h, h; z) = To(ki, k2; z)To(k3, k4; z)N^^\h, h, k 



4i ■ 



where 



N^'\h,h,h,h;z) = ([t;«(k3)0(k4)+t;(^)(k4)0(k3) + 5(k3,k4)] 



X Ro{z)[v^'\ki)(j){k2) + i;(^)(k2)0(ki) + 5(ki, ka)]) . 
After a significant amount of algebra we have the exphcit results for A^*^^^: 
ArW(ki, k2, k3, k4; z) = Nsubil2; 34) + Nr{12; 34) + NoiU; 34) , 
where we must symmetrize 

Ar^,Z5(12; 34) = NR,Dil2; 34) + iV^,,5(21; 34) + NR,Dil2; 43) + iV^,D(21; 43) 

with 

ju n9 S4^ = [ ^^^^ 

' ^ y (27r)'^ (27r)'^ (27r)'^ (27r)'^ 

X y(k3, k7, k8)\/(ki, k5, k6)2ro(k5, kg, k2)C'(24)C'(57)C'(68) 

where we use the notation C(68) = Cik^, kg), 

d'^h d'^h d% d'^ks 



{27rY {27tY {27tY {27rY 
X l^(k3, kr, ks)V{k^, k5, k6)2To(k5, kg, k2)2C'(27)C'(46)C'(58) 



and 



A^.«fe(12;34) 



d'^k. 



-4)ro(k5)\/(k5, ki, k2) V(-k5, k3, k4) 



{27lY 

X C'-l(k5)C'(ki)C'(k2)C'(k3)C'(k4) . 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 



(78) 



Put Eqs.(76), (77) and (78) into (74); in turn put EqJ74h into Eq.(72) and Eq.(72) back 



into Eq.(58) with n = 4 to obtain an explicit expression for C*-^-*(qi, q2; . 



Thus we have explicit expressions for C*^"^(qi, q2, -z) for n < 5. We use these results 
below. 
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V. EVALUATION OF MEMORY FUNCTION IN PERTURBATION THEORY 

A. Static part of Memory Function 

We want to determine the memory function K in a perturbation theory in powers of 
Di. We find that the static part of the memory function is of zeroth order in Di, while the 
dynamic part of the memory function begins at second order in Di. 

The static part of the memory function is determined by the equihbrium average 

J d^w if(^)(x, w)C'(w,y) = r('^)(x,y) = (50(y)^D<^50(x)). (79) 

In evaluating this static average it is very useful to use the identity: 

5B „ , . 5A 



{BD^A) = (3-' I rf'^xirf'^X2(^---r^(xi,X2)— -) (80) 

J ^0(Xi) ^0(X2) 

and we obtain 

r(^)(x,y)=^r'(r^(x,y)) . (81) 
It is easy to show, using Eq.([8]) for T^, that 

rW(x,y) = ■ V, [D5(x - y)] (82) 

where the average diffusion coefficient is given by 

D = (D(0)) = Do + ^i0o (83) 

where 0o = (0)- 



Taking the Fourier transform of Eq.(82) and multiplying by C ^{k) gives the static part 
of the memory function: 

K'^'\k) = t(3~^eDC-\k) = ik^Dx'\k) = iLo{k) . (84) 



Putting this result back into Eq.(14), dropping the dynamic part of the memory function and 
inverting the Laplace transform, gives the zeroth order approximation for the density-density 
time correlation function 

Co{k, t) = e-'='^^"'('=)*C'(fc) = e-^o('=)*C'(A;) (85) 

which agrees with the lowest-order result found previously. 
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B. Dynamic Part of Memory Function 

The dynamic part of the memory function for the dynamic structure factor is the sum of 
two pieces: 

r('^)(qi,q2;2;) = f(qi,q2;2:) + r,„b(qi,q2;z) (86) 
where the direct contribution is given h^p^ 

r(qi,q2;;.) = -{v\q2)Rizy (q,)) (87) 



and the subtraction part is given by Eq.(20) with 

iy(qi,ki;;.) = (0(ki)i?(^)t;'(qi)) (88) 

and 

v\q,) = tD'^<P{k,) . (89) 

We see that the dynamic part of the memory function vanishes at zeroth and first order in 
D,. 

We show here how to evaluate T^'^^ in perturbation theory up to fourth order. 

C. Direct contribution 

Focussing first on the direct contribution to the memory function we have 

d'^ki d'^k2 d'^h d'^ki 



r(qi,q2;^) = - J 



{27tY {2ttY {27rY {27rY 
X r(q2, kg, k4)V(qi, ki, k2)M(ki, k2, kg, k^; z) (90) 



where M is defined by Eq.(59). In perturbation theory 

d'^ki d'^k2 d'^ks d'^k^ 



f(")(qi, 



{27tY {2ttY {27tY {27rY 
\/(q2,k3,k4)l^(qi,ki,k2)M(")(ki,k2,k3,k4;^) . (91) 
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At second order M^'^^ is given by Eq.(61). Using Eq.(20), this reduces to 
f(2)(qi,q2;;,) = - 



d'^ki d'^k2 d'^h d'^h 



{27rY {27ry (27r)'^ {27rY 



V{q2, ks, k4)V^(qi, ki, k2) 



xTo(ki,k2;^)2C'(ki,k3)C'(k2,k4) . 

The third-order contribution vanishes, while the fourth-order contribution is given by 

d'^ki d'^k2 d'^ks d'^k^ 



(92) 



f(^)(qi,q2;^) = -/ 



{2txY (27r)'^ (27r)'^ (27r)' 



X V{ci2, ks, k4)To(k3, k4; z)\^(qi, ki, k2)To(ki, k2; z) 



X 



iV(4)(ki,k2,k3,k 



4; Z) 



(93) 



is given by Eq.rt74|. This contribution divides naturally into three pieces. One piece, 



ritb' when added to T^^"!,, vanishes. We have then that the fourth-order contribution is given 
by 



r(^)(qi, q2; z) = f g)(qi, q2; z) + f ^'\qi, q2; z 



(4)/ 



(94) 



where 



r2^(qi,q2;2;) = - 



d'^ki d'^k2 d'^ks d'^k^ d'^k^ d'^ke d% d'^ks 
{2nY {2nY {2ny {27tY {27tY {27iY {27iY {27iY 



X V{q2, ks, k4)To(k3, k4; z)V{qi, ki, k2)To(ki, k2; z) 



X 4\/(k3, kj, kg) V(ki, ks, k6)2To(k5, kg, k2)C'(24)C'(57)C'(68) 



(95) 



and 



rDHqi,q2;2;) 



d'^ki d'^k2 d'^ks d'^k^ d'^k^ d'^k^ d% d'^kg 
{27iY {27iY {27iY {27iY {2nY {27iY {2nY {27iY 



X V{q2, ks, k4)To(ks, k4; z)V{qi, ki, k2)To(ki, k2; z) 



X 4V{k,, ks, k6)l^(ks, k7, k8)2To(k5, kg, k2)2C'(27)C'(46)C'(58) 



(96) 
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VI. BARE PERTURBATION THEORY AT SECOND ORDER 

A. General Form 

Now that we have the perturbation theory results we need to see the physical conse- 
quences. We begin with bare perturbation theory at second order. For a general static 



structure factor, T^'^\q, z) is given by Eq.(92). After integrating over the 5-functions we 
have 

where 

y(q,k) = '-D,q ■ [kx-\k) + (q - k)x-^(q - k)] . (98) 



B. Structureless approximation 

In the structureless approximation we assume that the static susceptibility is independent 
of wavenumber 

X-^(k) = r 

and introduce a large wavenumber cutoff A. This approximation (model) is appealing for 
two reasons. First, in this case, the vertex simplifies to the form 

V{q,k)='^D,rq' . (99) 

Second this model corresponds to a coarse-grained system where one has integrated out 
short-distance degrees of freedom including the first peak in the static structure factor. 



Eq.(97) then becomes 



r(^)(,..)^lz.;,V-/(g, ^,^.,^^,,.V(^_,)., . (100) 

Letting k = | + p in the integral gives 

1 rl'^n 1 



{27rY z + iDr[2p^ + q"^ /2] 
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Doing the angular integral 

r'^'fc.) - ^^?''r^ir./;/-Mp^^^5^pL^ (101) 

where 

/rid I. 
. (102) 

At this point we move to dimensionless variables. If we set p — Ax in the integral then 

r«te..) = l^'^^nA^ (103) 

where Q — q/A and we introduce the time r = Then we have 

T^'\q,z) = -i^Q^No{n) (104) 

T 

where 

No(n) = KJai^) (105) 

and 
with 

n = -i^r + . (107) 

The dimensionless coupling g is given by 

^^^'^'s (108) 



2 V L> 

where 



S^CA'^^ {{Spf) . (109) 

To see that g is dimensionless, note that Di has dimensions of £>/0o where 0o is the equi- 
librium particle density which has dimensions of A*^. Finally the fourier transform of the 
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static structure factor has dimensions C ~ L'^Pq ~ A'^. Because of the factor in Eq.(104) 
we see that there is no second-order contribution to the diffusion coefficient. 



To go further we must evaluate the dimensionless integrals in Eq. ( 106 ) . In two dimensions 
we find 



In three dimensions we have the explicit result 



1 



-tan 



-1 



In the small q and z limit we have for general d: 

IJQ) = - ^'^"'^dx = , } . 
^ ' 2 Jo 2{d-2) 



which is well defined for ci > 2 and 



r(2)(g,o) = -zVc' 



T 2{d-2) 

The kinetic equation in bare second-order perturbation theory is given by 

'z + zLo{q) + K'^'\q,z)]C{q,z) = '^ . 
L J r 

This can be written in dimensionless form 



[u + iQ^D{Q,u)]C{Q,u) 



where i' = zr and the damping is given by the real part of 



The dynamic structure factor is given by 



S{Q,v) = -2nlm 



'/3-W 1 
r ly' + iv" 



2np-W u" 

r (iy')2 + (zy")2 



where 



(110) 



;iii) 



(112) 



(113) 



(114) 



(115) 



(116) 



u" = Q'{l-gQ'KJ') 



;il7) 
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FIG. 1: Plot of the dynamic structure factor versus wavenumber for zero frequency for different 
values of the coupling: g=50, 60, 70, 80 and 85 from bottom to top. 



jy' = u{l + Q'gKdJ) . 
For d = 3 we have the simple integrals 



/' = / x'^dx 
Jo 



Z/2 + X? 



Jo 



z/2 + xf 



(118) 

(119) 
(120) 



where xi = ^ + 2x^. 

We plot the dynamic structure factor in Figjljfor d = 3. The conservation law dominates 
the structure for small wavenumbers. However for large wavenumbers one sees the devel- 
opment of an instability. In this approximation the instability comes from short distances, 
Q = 1, and low frequencies as seen in Figjl] 

In Fig|2]we plot the dynamic structure factor at Q = 1 versus frequency. The instability is 
manifest at low frequencies. 
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FIG. 2: Plot of dynamic structure factor versus frequency for Q = 1 for different values of the 
coupling: g=65, 70, 75, 80, and 85 from bottom to top at u = 0. 

The instabihty in the problem sets in for the couphng g where ReD first becomes negative. 
In three dimensions the critical g is given by 

1- g*K^r{Q = l,u = 0)=0 . (121) 

This has a solution 

27r2 

9* = —. — ^ ■ (122) 



We have from Eq. ( 70 ) 



I\Q = l,v = Q) = ]^[l-han-\2)] . (123) 



Given that tan-^{2) = 1.107. .. we find g* = 88.42 .... 
In two dimensions 



X 1^ xdx 1 /2 + Q2\ 
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and the damping is given by 

:)2 = 1 - r7(02_ 

in \ Q 



D{Q\0) = l-gQ'j-ln('^-±f-] . (125) 



which leads to the instabihty couphng 



C. Gaussian Structure in Hydrodynamical limit 

How strongly does our result depend on the cut off A. To see this consider the case where 
xiq) falls smoothly to zero for large q. It is simplest to look at the problem in the small q 
and z regime. In this case, to leading order in q, the vertex takes the form 

F(q,k) = '-D^x-\k)q,q, [5,, - hk,a{k)] (127) 

where sums over i and j are implied, 



and from Eq.(97) for small q: 



r(2)(g, z) = -^qiqjqkqm / [Stj - hkja^k)] [5km - kkkma{k)] . (129) 

A practical choice for the static susceptibility is given by 

X{k) = Xoe-5(^^)' (130) 
where i is the characteristic length and 

a = ^ = -e^ . (131) 



In the z = limit, after considerable algebra, Eq.(129) reduces to 

T^'\q,0) = -z~gQ'f~'C{0)j, (132) 

where 

_ K,^id-mr (^] + ^^d' - n33) 
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where 

^ (134) 

and Q = qi. We see that the resuhs for r^^)(g, 0) are very similar for the two different 
choices for static structure factor if we make the correspondence A ^ 1/i and r Xo^- 



VII. BARE PERTURBATION THEORY AT FOURTH ORDER 

Here we look at the reduction of the two-loop contributions to the memory-function in 
more detail. There are two contributions. 



A. General Reduction of T 



(4) 
D 



We have from Eq.(96), 



d'^ki d'^k2 d'^h d'^h d'^h d'^h d% d'^h 
(2tiY {2tiY {2txY {2tiY {2txY {2tiY {2txY {2t:Y 



X V{ci2, ka, k4)To(k3, k4; z)l^(qi, ki, k2)To(ki, ks; z) 



X 4V^(ki, ks, k6)V^(k3, k7, k8)2To(k5, kg, k2)2C'(27)C'(46)C'(58) . (136) 

First do the integrations over the 5-functions associated with the C, then over those associ- 
ated with the cubic vertices. This leads to the result 




X To(k3, -qi - ka; 2;)qi ■ A(ki, qi - ki)To(ki, qi - ki; z) 
X ka ■ A(qi - ki, ki - qi - ka)ki ■ A(ki - qi - ka, ka + qi) 



xTo(ki-qi-ka,ka + qi,qi-ki)C'(qi-ki)C'(-qi-k3)C'(ki-qi-ka) . (137) 
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B. General Reduction of T 



(4) 
R 



We have from Eq.(95)) 

rg^(qi,q.;.) = - 



d% f d'^k2 f d'^h f d'^h f d'^h f d'^k^ f d% r d'^k^ 



{2ttY J {2ttY J {2nY J {2ttY J {2ttY J {2ttY J {2tiY J {2tiY 
X V(q2, ks, k4)To(k3, k4; z)V{c[x, ki, k2)To(ki, k2; z) 



X 4r(k3, kr, kg) V(ki, ks, k6)2ro(k5, ke, k2)C'(24)C'(57)C'(68) 
Doing the integrations over the internal (5-functions leads to the result 



(138) 



rg)(qi,q2;^) = - 2l^^(2vr)'^5(qi + q2) ^ ^^^^^^ 



1 



d'^ki 



qi • A(ki,qi - ki' 



X ro2(ki, qi - ki; z)C{qx - ki)7(ki, z + zLo(qi - ki)) 



(139) 



where the insertion 7 is defined 



7(ki,2;) 



d'^k^ 
{2tiY 



ki-A(k5,ki-k5) To(k5,ki-k5,qi-ki)C'(k5)C'(ki-k5) .(140) 



C. f In the structureless approximation 



In the structureless approximation the interaction vertices simplify significantly and 



Eq.(137) becomes 



ri,^)(qi,q2;z) = (27r)''<5 (qi + q2) r^^^(qi, ;^ 



(4)/ 



(141) 



where 



{2ny J {27rY 

X To(ki, qi - ki; z)To{ki - qi - ks, ks + qi, qi - ki) . 

As a first check on this result let us look at the small qi and z limit where we find 

,C 



where 



J, 



4{d - 1) Jo 



y"'-'dy 



2y+l 2 + \y 



(142) 



(143) 
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D. In the structureless approximation 



Similarly we evaluate given by Eq.(139) in the structureless approximation and 
obtain in the long time and distance regime 

f'i\q^,ci,;z) = (2vr)^5(qi + q2)fg)(qi;^) (144) 

and 



T 2 



where Jd is given by Eq.(143). 



(145) 



E. Summary of Bare Perturbation Theory Results 



Combining the small q and z limits for terms up to fourth order we have 



T^''\q,^) = -igQ'-Kdl 

T 



where 



7 



1 3 ~ 

+ -gJdKd 



2(rf-2) 2 

One can interpret this in terms of an effective coupling 

geff = g{l + 'i{d-2)gJdKd) . 
For perturbation theory to make sense we require that the coefficient 

Cd = ?>[d-2]JdKd 



be small. In three dimensions C3 = |^ where 



1 /"l 

J?, = - I dxx^ 



A Jo 



2 2 

+ 



4a;2 + 3 4 + 3x'^ 



0.0622. . 



and C, = 0.00946. . .. 



VIII. SELF-CONSISTENT PERTURBATION THEORY 



(146) 



(147) 



(148) 



(149) 



A key ingredient of MCT is that it is a self-consistent theory where the memory function 
is a function of the full correlation function. Here we show how this is arranged through 
two-loop order in our development here. Elsewhere we discuss how this is naturally carried 
out in the MSR formulation. 
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A. Second Order Theory 



One wants to replace bare correlation functions by renormalized correlation functions. 



We begin with the bare second-order memory function given by Eq.(92) 

d'^h d'^k2 d'^ks d'^ki 



{27:)^ {2'kY {2ttY (27r)' 



X 



X V^(qi, ki, ks) V(q2, ks, k4)ro(ki, ks; z) 



^(ki, k3)C'(k2, k4) + ^(ki, k4)C'(k2, k3 



We first write this in terms of the bare two-point correlation functions. We have 



(150) 



ro(ki, ks; Z) ^(ki, k3)C'(k2, k4) + ^(ki, k4)C'(k2, ks) 



-% I (ite'^^e-^oe^l^^e-^"^'^)* ^(ki, k3)C'(k2, k4) + ^(ki, k4)C'(k2, ks) 



/ dte'^' [Co(ki, kg; t)Co(k2, k4; t) + Co{k,, k^; t)Co(k2, kg; t)] . (151) 
Jo 



At this order we can replace Cq C, and the last equation is replaced by 



-t I dte'^' [C(ki, kg; t)C(k2, k4; t) + C(ki, k4; t)C(k2, kg; t)] 
/o 



and we have for the memory function at second order 

d'^ki d'^k2 d'^ks d'^k^ 
(27rp(27rp(27rp(27r)^ 



rgHqi,q2;^) = -/ 



\/(qi,ki,k2)F(q2,k3,k4) 



(152) 



X 



-t I dte'^' [C(ki, ka; t)C(k2, k4; t) + C(ki, k4; t)C(k2, kg; t)] 





(153) 



This form will generate contributions at 4th order in perturbation theory. Generating the 
4t/i-order contribution from this result, requires generating the second-order contribution to 



the correlation function. Iterating Eq. ( 14 ) 



C(qi, q2; z) = To(qi; z)C{qi,q2) - To(qi; z) 



d% 
(2^ 



K(2)(q^,ki)C(ki,q2)ro(q2; 



To(qi; 2;)C(qi, q2 - To(qi; 2;)r(^)(qi, q2)To(q2; z) 



Taking the inverse Laplace transform to go to the time domain gives 

C(qi,q2;t) = e-^«(^)*C'(qi,q2) 



JO Jo 



We then substitute this result into Eq.(152) and keep terms of Ath order. We find 



rt 

X 



(qi, q^; t) = j (^^(qi' ki, k,)l^(q2, ka, k4) 

POO 

X (2z) / rfte^"*[2Co(k2,k4;t) 

JO 

X rtise-^«('^i)(*~^) rrfrr(2)(ki,k3,s-r)e-^°(^^)^ 

JO JO 

/■CO 

X (2i) / cite*"*e-^°(''^)*2Co(k2,k4) 

JO 

r cise-^»('^^)(*-'^) r rfrV^(ki, k5k6)V^(k3, kykg) 
JO JO 

X e-(^«('^5)+^«('^«))(^-^)2C'o(k5, k7)C'o(k6, k8)e-^"(^^)") 
= I (2^ (2^ (2^ (2^^^^^' k4)2Co(k„ k4) 

X {2i) / cite*"* / rfs e-(^o('^i)+^°(''2))(*-^V(ki,k5,k6)\^(k3,k7,k8) 

JO JO 

X (ire-(^o('^5)+^°(''«)+^«('^2»("-")2C'o(k5, k7)C'o(k6, kg) 

JO 

^ g-(Lo(k3)+Lo(k4))r) 

/ 4^4^4^i^^('li''^i''^2)^('l2,k3,k4)Co(k^ 



X y(ki, kske) V(k3, k7k8)ro(k5, kg, ks; z)2C'o(k5, k7)C'o(k6, k; 



which agrees with rj^''(qi, q2, 2;) given by Eq.(95). So r^"* is generated by expanding F 



.(4) 
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B. Two Loop self-consistent theory 

We want to replace the Ath order bare contribution with a self- consistent form which 
depends on the full correlation functions. We begin with the bare contribution 

d'^ki d'^k2 d'^ks d'^k^ d'^k^ d'^k^ d% d'^ks 



fg^(qi,q2;^) = -16 / 



{2nY {2nY {2nY {2nY {2nY {2nY {2nY {2ny 
X y(qi, ki, k2)ro(ki, k2; z)V{c[2, ks, k4)To(k3, k^; z) 

X V{k,, kg, kg) V^(k3, ky, ks)To{k,, ke, k2)C'(27)C'(46)C'(58) 

_ r d'^ki d'^k2 d'^ks d'^k^ d'^k^ d'^ke d% d'^ks 
J {2nY {2nY {2nY {2nY {2nY {2nY {2nY {2nY 

X V{q,, k^, k2)V{q2, ks, k,)V{k^, k,, k^Mk^, k^, kg) 

X C{27)C{46)C{58) 



X 



[z + iLo(l) + iLo(2)] [z + iLo(3) + iLo(4)] [z + iLo{2) + iLo(4) + iLo(5)] 

We can then use the following result based on the pole structure of the zeroth order 
correlation function: 

rfcui C„(/,-i.^'i) C'(A-i) 



/ 



27r 2; - u;i + iLo(2) 2; + iLo(l) + iLo(2) 
Using essentially this result five times we find 

d'^ki d'^k2 d'^ki d'^k^ d'^k^ 



fg^(qi,q2;^) = -16 / 



{2t:Y {2t:Y {^ttY (27r)'^ (27r)<^ 
X V{q2, k3, k4)V^(qi, k^, ks) V(ki, k,, -k4)y(k3, -ks, -k,) 

xC-\k,)C-\k^) ^ ^ ^ 



[z-Ui- U2] [z-ljJ2,- OJ^ [z -U2-OO4- 005] 
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To obtain the self-consistent generalization to this order we replace 

Co^C . (154) 

We then have 

d% cfik2 cfih d'^h cfih 



rg^(qi,q2;^) = -16 / 



{2t:Y {2t:Y {2-nY {2t:Y {2t:Y 
X V^(q2, ks, k4)y(qi, \i^M)C-\h)V ks, -U)C-\h)V -ks, -K) 

X j ^^CR{ki,z - uj2)C{k2,uj2) j ^^CR{k2„z -uJi)C{ki,uJi)CR{kf„z -uj2-uji) 
where the retarded correlation functions are defined by 

C«(*,.-)=/^^i^ (165) 

If we define a vertex 

t/(qi,ki,k2) =C'-i(?i)y(qi,ki,k2) 

then 

d'^ki d'^k2 d'^k-s d^^k^ d'^k^ 



rgHqi,q2;^) = -16C'(qi)C(q2) J 



{27rY {2nY {2nY {2nY {2nY 
X t^(q2, kg, k4)V^(qi, ki, k2)l^(ki, ks, -k4)V^(k3, -k2, -ks) 

X J ^CR{ki, z - uj2)C{k2, UJ2) j ^CR{k3, z - u;4)C(/c4, a;4)Cij(A;5, z-u;2- (^4) 

The most useful form for our purposes is 

p(4). . . f d'^ki d'^k2 d'^kz d'^ki d% 

ry(q:,q2;.) = -16C(q,)C(q2) ) (2^(2^(2^(3^(2^ 

X y(q2, ka, k4)y(qi, ki, ks)^ (ki, ks, -k4)y(k3, -^2, -ks) 

X / ^C(fc„c.O / '^C{k2,.2) J '-^Ciks^us) I ^C(A;4,.4) / ^C(fcs,a;5) 

11 1 

X 



[z-Ui- U2] [z-(jJ2,- OJ^ [z -U2-OU4- ^5] 
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We need to invert the Laplace transform and obtain this contribution in the time regime. 
The key result we need is 



/ 



1 1 1 

27ri [z — LUi — LJ2] [z — U!3 — UJ4] [z — UJ2 — 1-^4: — (^s] 



Jo Jo 

which is a product of convolutions. One can the do the frequency integrals easily to take 
one fully to the time regime: 

X V(q2, k3, k4)y (qi, ki, ks)^ (ki, ks, -k4)y(k3, -ks, -kg) 

X /* dh r dt2C{ki,t - h)C{k2, t - t2)C{k3, t2)C{k4, h)C{k^, h - is) ■ 
Jo Jo 

In the structureless approximation this reduces to 

rgnqi,q2;t) = {2nr6{{ci, + ci2)T^^\qut) (156) 
and, in terms of dimensionless variables, 

rS^qi; t) = 4^CQ^ / jl dt, dt2f{K,, T - T^Q - K,,T - T,) 

xf{K3,T2)f{Q + K3,T^)f{-Q + K,-Ks,T^-T2) . (157) 

We show elsewhere that this same structure is found in the time regime at two-loop order 
using the MSR formulation. 
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C. Self-Consistent Kinetic Equation at Second Order 

The second order memory function in terms of full correlation functions is given by 

POO 

Jo 

In the structureless approximation this reduces to 

T^i\q,t) = K('^'\Q,t)C (158) 



and 



K^'''\Q,t) = Q''\No{K,t) (159) 



r2 



where 



No{Q,T)=g J ^J{K,T)fiQ-K,T) 



where 



(2vr) 



f{Q,T)=C{Q,T)/CiQ) (160) 



and we have introduced the same dimensionless variables as in the case of bare perturbation 
theory. The kinetic equation reduces to 



df{Q,T) 
dT 



Q'fiQ, T) + r dSNoiQ, T - S)f{Q, S) . (161) 
Jo 



We first look at the solution for f{Q,T) numerically. We begin, in Figjsj with small g and 
find, as in bare perturbation theory, near exponential decay with time for fixed wavenumber. 
Another way of characterizing the data is in terms of a running relaxation time 

r(Q,T)= r dtf{Q,t) . (162) 
Jo 

This quantity is shown in Figjl] for a series of short times and coupling g = 10. Notice 
that t{Q,T) has approached t{Q) = t{Q, 00) for the largest times shown for Q > 0.75 
but not for smaller Q. If we increase g to 60.0 we see substantial slowing down as shown 
in FigjS] These trends are continued as we move to g = 90.0 as shown in Figj6] with the 
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Q 



FIG. 3: Plot of normalized intermediate structure factor versus wavenumber for different times 
T = 0.1,5.0, 10.0, 15.0, and 20.0 from top to bottom. The coupling g = 10.0 is relatively weak. 

new feature, shown more clearly in FigjT], that we see the development of a weak peak near 
Q ~ 0.80 which saturates at T = 200.0. Finally, with g = 100.0 the system, shown in FigjSj 
is rendered unstable and a peak at Q = 0.88 grows rapidly with time in the structure factor. 
We show in Fig|9]that this large wavenumber peak can be fit to a gaussian form 

/p(Q,T) = Ae-^(«^««)' (163) 

where A is the peak amplitude, B is a new growing length squared in the problem and Qq 
is an ordering wavenumber. 

The central Q = peak can also be fit to a gaussian with Qq = and A = 1 and another 
growing length squared Bq. We find Bq ^ B and both grow as a power law with time. 

We see that as g goes from 90 to 100 that the system goes from stable to unstable. How 
do we find the transition value (?*? A good way of determining g* is to work in the unstable 
phase where the structure peak amplitude A has a minimum at time T*, as a function of 
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FIG. 4: Plot of T(Q,r) versus wavenumber for different times T= 0.1, 5.0, 10.0, 15.0, and 20.0 
from bottom to top. The coupling g = 10.0. 



T as shown, for example, in Fig, 10 The closer g is to g* the longer T*. If one plots, as in 



Fig, 11, T* versus g and fits T* to a power-law diverging aX g = g* one obtains an accurate 



estimate of g* . A good fit to the data shown in Fig, 11 gives an estimate g* = 93.0047558. If 
we plot the value of the amplitude minimum, Amm, for g > g*, we see in Fig 12 a.sg^g* 
from above, Amin appears to go to zero. 

Next we look at the behavior of the intermediate structure factor peak for g = g*. In 



Figs 13 and 14 we plot the structural peak for intermediate times and show the gaussian 



fits. The fit parameters, Qo, B , and A, are shown as functions of time in Figs|T5| [16] and 
17 Qq orders rapidly, while B can be fit to a simple power- law form. For g < g*, A can 
reasonably be fit using 

„-ET 



A = A, 



(T + to) 



(164) 



Such fits are shown in Figs, 17 and 18 The fit in Fig, 17 is over a very long time scale and is 
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FIG. 5: Plot of t{Q,T) versus wavenumber for different times T= 10.0, 20.0, 30.0, 40.0 and 50.0 
from bottom to top. The coupling g = 60.0. 



breaking down for the longest time. Carrying out fits using Eq.(164) for a range of values 



of the coupling constant, we find the parameter Aq shown as a function of g in Fig, 19 The 



parameters E and a are more interesting. E, as a function of g, shown in Fig, 20, can be fit 
to the form 



XE 



E = AE{g*-g) 



with fitted parameters Ae = 0.0137, xe = 0.768 and g* 



(165) 

93.0067. It is clear physically 



16 



that the critical point corresponds to E 0. The exponent a is plotted in Figj2T 

The new kinetic length in the problem is ^/B. As shown in Fig 
B can be fit to the form 

B = BnT'"' 



versus g. 



for the critical coupling, 



(166) 



The exponent xb is shown as a function of g in Fig, 22 Away from the critical point B 



grows nearly linearly with time but it crosses over to 1/3 as g —>■ g* . 
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FIG. 6: Plot of r(Q, T) versus wavenumber for different times T= 10.0, 50.0, 100.0, 150.0 and 
200.0 from bottom to top. The coupling g = 90.0. 



One can return to the relaxation time defined by Eq.(162) with t{Q) = t{Q, oo). In 
Fig, 23 we plot the the maximum {Q > 0.5) of t{Q) versus coupling g. We see a very strong 
dependence on g as it approaches g*. 
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FIG. 7: Plot of r(Q, T) versus wavenumber for different times T = 10.0, 50.0, 100.0, 150.0 and 200.0 
from bottom to top. The coupling g = 90.0. 



One can estimate t{Qq) using Eq.(164) for the peak contribution to the structure factor, 

-Et 



poo 

r{Qo) = dtAo 
Jo 



e 



{t + toY 



POO 

= AoE"-^ / dy- — . (167) 

Jo \y + toEY ^ > 

Assuming as ^ 0, that IqE goes to zero, we have 

r(go) = ^o^""'r(l - a) . (168) 

The results E ^ [g* — g)^^^ and t{Qo) ^ {g* — g)^^^^ imphes that a = 5/9 which is 



compatible with the value for a given by the fit in FigjTT 

We conclude that in this model there is the development of unanticipated structure in 
the intermediate structure factor for g near g*. In the stable regime the associated peak 
narrows thus giving a growing lengthV^ in the problem. In the next section we show how 
some approximate analytical progress can be made looking at this peak formation. 
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FIG. 8: Plot of normalized intermediate structure factor versus high wavenumber for different 
times T = 5.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0 and 40.0 from bottom to top at Q = 1. The coupling 
g = 100.0 is above the critical value. 

IX. PEAK AMPLITUDE EQUATION 

A. Mapping onto Amplitude Dynamics 

In order to make analytical progress on our one-loop direct theory we assume that our 
long-time solution is of the form 



f{Q,T) = UQ,T) + f,{Q,T) 



(169) 



where 



MQ,T) = e 



-BoQ^ 



(170) 



and fp is given by Eq.(163). B and Bq are large and A grows or decreases with time 



depending on whether we are in the stable or unstable phase. Qq is a fixed wavenumber 
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FIG. 9: Plot of the normalized intermediate structure factor versus wavenumber for time T = 40.0. 
The coupling g = 100.0 is above the critical value. Shown also is a fit to a gaussian of the form 



given by Eq.(163). 



characterizing the position of the peak. Let us define 



TT 



Bo{T) 



and 



(171) 



TT 



B{T) 



Eq.(169) is then of the form 



(172) 



f{Q, T) = 0o(r) (Q, T) + 0(r) As(Q - Qo, T) 



(173) 



where 



(174) 
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FIG. 10: Typical behavior of peak amplitude A versus time for coupling g > g* . This example is 
for g = 93.005. 



and 



lim As (Q,T) =5(Q) . (175) 

T — >oo 



We then substitute the assumed solution, Eq.(169), into the kinetic equation 

^^li9jIl = _Q^f^Q^T) + Q''£dSNiQ,T-S)f{Q,S) . (176) 

and look for self-consistency. A key assumption is that the length squared B is arbitrarily 
large. We work here in three dimensions. 

The peak near Q = is simpler to treat due to the explicit Q dependence in the kinetic 
equation. For small enough Q one can drop the interaction term and one has 



dfo{Q,T) 



dT 

with the solution 



Q'fo{Q,T) (177) 



/o(g,T) = e^«'^ (178) 
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FIG. 11: Plot of the time T* where the peak ampHtude A hits a minimum before going unstable 
versus g. There is an excellent fit to T* = 435445/(5 - c/*)0-6389 ^^^j^ g* ^ 93.0047558. 

which gives Bq = T. 

Next we focus on the peak near Qq. We can write 



dfp{Q,T) 



dT 

and we need to evaluate the memory kernel 



QlfpiQ, T) + Ql r dSNpiQo, T - S)UiQ, S) (179) 

t/ 



(2 

for Q = Qo- Assuming the 5-function form 



N,{Q, T) = gJ j^.fpiK, T)/,(Q - K, T) (180) 



f,{Q,T) = cP{T)6{Q-Qo) , (181) 



one can do the K itegration in Eq.(180) with the result 



N,{Qo,T) = ^^<P\T) . (182) 
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FIG. 12: Plot of Amin, the peak amplitude minimum before going unstable, as a function of g. 
We are left with the kinetic equation valid near Q = Qo 

^^'f^ = -QlfpiQ, T) + Ql dSN,{Q,, T - S)f,{Q, S) . (183) 

CanceUng a common factor of the 5-function, gives the equation for the peak amplitude 
equation 

0(T) = -Ql(j){T) + GQl r dS<p\T - S)<P{S) (184) 

JO 

where 

G=|| . (185) 



Changing the scaling of time to t = Q^T we obtain 

0(t) = -0(t) + G f ds(f)\t - sWs) . (186) 
Jo 

If one replaces 0(s) by 0(s) inside the integral this equation of motion reduces to 
Leutheussar's equation!^. 
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FIG. 13: Plot of the normalized intermediate structure factor for high wavenumbers at the critical 
coupling g = g* for T=500, 1000, 1500, and 2000. Shown also is the gaussian fit to the structural 
peak. 



We will assume that Eq.(186) can be solved as an initial-value problem with (f)(0) = 1. 



We find that self-consistently that we have been able to replace the kinetic equation, 



Eq.(161), with the amplitude equation. Eq.(186). We now show that the solutions to 



Eq.(186) shows the same phase structure as found in the numerical solution of Eq.(161). 



B. Power-Law Solution 



In terms of Laplace transforms the equation of motion for the amplitude satisfies 



[z + i + N{z)]4>{z) = 1 



;i87) 



where 



[z) = -i dt e'''(j){t) 
Jo 



;i88) 
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FIG. 14: Plot of the normalized intermediate structure factor for high wavenumbers at the critical 
coupling g = g* for T = 2000,4000,6000,8000, and 10000. Shown also is the gaussian fit to the 
structural peak. 

and 

roc 

N{z) = -iG / dt e'^'^^t) . (189) 
Jo 

We want to show in the long-time limit and near the critical coupling G* there is a power-law 



solution to Eq.(186) of the form 



We want to determine the exponent a. 

The first step is to work out the Laplace transforms for (p and for the trial solution. 
We have 

An 

(h(z) = -i dt e'^' ° 
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FIG. 15: Plot of the ordering wavenumber (X 100) versus time at the critical coupling g = g*. 

^ r ^e^-(-*o) _ (191) 

J to X 

Let y — xz m. the integral to obtain 

(j){z) = -aoe"^*°'^""V(a) (192) 

where the integral reduces to 

(^\0i) = / — e^ = cro-— ; :- + ... (193) 



and 

Next look at the memory kernel given by 



oo 4^ 



N(z) = -i / dt e'^'G ° 

^0 
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FIG. 16: Plot of the peak width versus time at the critical coupling g = g* ■ Fit is shown to a 
power law form B = 16.21(r + 121.2)0-293. 



"^"*o[A^(0) + AiV(2)] 



(195) 



where 



A^(0) = -lAlG 



dx 



to X 



2a 



-iAlG 1 



2a - 1 t^Q 



2a- 1 



and we have assumed that a > 1/2. In Eq.( |195[ ) 

AN{z) 



-iAlG 



dx 



to 



2a I 



1] 



(196) 



(197) 
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FIG. 17: Plot of the peak amplitude versus time at the critical coupling g = g* ■ Fit shown is to 



the form given be Eq.(164) with Aq = 0.672, E = 4.70 x 10"*', to = 9.00 and a = 0.55. 



To lowest order in z the power-law solution corresponds to the cancellation of terms in 
the kinetic equation 

+ i + Ar(0) = (198) 

which gives the result 

^ -1 (199) 



2a - 1 tl""-^ 



which depends explicitly on the time cutoff. We turn to the next order term in the small z 
expansion of i + N{z) which is given by 



AN{z) = -^Alz^--'G ^^^-1] 

Jtoz y 

= -iAlz^^'^GJ{a) (200) 
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FIG. 18: Plot of the peak amplitude versus times at coupling g = 93.004720. Fit is shown to the 



form given by Eq.(164|. The fit is indistinguishable from the numerical data. 



where we have the remaining integral 



oo 



^[e^^-l] = Jo(«)-^%f^ + ... (201) 
Jtoz 2/^" 2(1 - a) 



where 



'^o(«)= r%[^''-^\ ■ (202) 
JO y 

Assuming 0.5 < a < 1 the integrals ctq and Jq can be evaluated. The kinetic equation 
then takes the form 

[z + i + e-^"*° [iV(0) - iAlz^^-^GJ{a)] [-aoe-^*«"z"- V(a)] = 1 . (203) 
Since z^"-~^ <^ z this reduces to 

[-iAlz^''~^GJo{a)][-iAQe-'''>'z^-^ao{a)] = 1 (204) 
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FIG. 19: Plot of the parameter Aq as a function of g. Fit is given by 0.194/(93.0048 - 5)°-^°^. 
which requires 

2a - 1 + a - 1 = (205) 
or q; = 2/3 which is a self consistent value. We are left with the equation 

- AlGJo{2/S)ao{2/3) = 1 (206) 
It is left to appendix C to show that 

(70 (2/3) Jo (2/3) = -2nVS ■ (207) 
We then have the constraint on the solution 

27iV3AlG = 1 . (208) 
Notice that the results for a and Aq do not depend on the short-time cutoff. 
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FIG. 20: Plot of the parameter E as a function of g. The fit is to the form given by Eq.(268) with 
Ae = 0.0137, XB = 0.768, g* = 93.0067. 



C. Numerical analysis: One-loop Case. 



We can numerically solve the amplitude equation, Eq.(186). We first determine the times, 



t*, when in unstable runs, hits its minimum versus G. A power-law fit assuming t* goes 
to infinity as G goes to G* gives an estimate G* = 0.79992 .... 

In the stable regime, G < G*, the peak amplitude decay can be fit to the form 



(f){t) = Ao 



{t + to 



(209) 



In Figs,24| 25, and 26 we plot the fit parameters Aq, E, and a as functions of G. 

Numerically we find at the critical point, where E = 0, Aq = 0.51, to = 0.185, and the 



exponent a = 0.668. The analytic result, Eq.(208), with G* = 0.79992 gives Aq = 0.508. The 



analytic results agree with the numerical results. We can also compute the relaxation time r 
as a function of G. We obtain a very good fit to the data with the form: r = tq/{G* — G)^^! 
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FIG. 21: Plot of the parameter a in Eq.(164), as a function of g 



with To = 1.387, G* = 0.799925 and x^, = 0.2197. For G = G* we find r(t) = where 
ri = 2.156 and x^^ = 0.2718. 



D. Two-Loop Amplitude Contribution 

We now work out the results of the projection onto the structural peak solution at two- 
loop order. We begin with the two-loop expression for the dynamic part of the memory 
function in the structureless approximation in terms of the dimensionless parameters intro- 
duced at one-loop order: 



X f{K,,T - T,)f{Q -K^,T- T,)f{K,, T2)f{Q + K^, T^)f{-Q + K^- K^, T, - TsplO) 



We again assume a trial solution given by Eq.(181) where the unstable wavenumber Qq is 



time independent. We restrict the analysis here to three dimensions. We need the memory 
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FIG. 22: Plot of the parameter xb in Eq.(268), as a function of g. 
function evaluated at Q = Qo and 



ir('^'^)(Qo,T)=4^Q^J(Qo,T) (211) 

where 

J{Qo,T) = dT, dT,<p{T - TMT - T,)<P{TMT,)<P{T, - T,) J ^(2^^i^l 

X S{K, - Qo)S{\Q - - Qo)SiK, - Qo)Si\Q + K,\ - Qo)Si\ - Q + - K^l - Qo)(212) 

One can then do the integrations over Ki, K^, ui = Ki ■ Q, u-^ = ■ Q and the azimuthal 
angles with the result: 

J(Qo,T) = ^ £dT, dT^(^{T - ri)0(T2)0(r - T^)cl>{T{)cl>{T, - T,) . (213) 

The contribution to the memory kernel at two-loop order is given by 

ir('^'^)(Qo,T) =4^^^ j^dT, J^'dT2<P{T~TMTMT-T,)<P{T^)<P{T,-n)i^^^^ 
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FIG. 23: Plot of t{Qo) as a function of g. Fit is to form To/{g* — gY'^ with values tq = 8.96, 
g* = 93.0048 and Xr = 0.344. 



The two-loop peak-amplitude model is given by 

d 
di 



(t>{t) = -QUit) + f dsN{t - (215) 



wiith the memory kernel 

H ft 



N{t) = GQt(p\t) + GiQl j dhT dt2(p{t - h)(P{ti)(P{t - t2)4>{t2)4>{ti - t2) (216) 

JQ Jo 

and we have the couplings 



We then rescale times t — Q^T and have 

d 

dT 



G = (217) 



G..4/%^^^ . (218) 



0(r) = -0(T) + r dsNiT - SUiS) (219) 
Jo 
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FIG. 24: Plot of ^o, defined by Eq.([209j), as a function of G. 

and 

N{T) = Gct>\T) + Gi r dT, r dT,<p{T - T,)(f){T,)<p{T - T2)0(T2)0(Ti - T^) . (220) 

Jo Jo 

This model can be solved numerically, first we look at the power-law solution at two-loop 
order. 

E. Power-Law Solution at Two-Loop Order 

We insert the trial solution (changing notation from T to t) 



into the two-loop contribution in Eq.( 220[ ) with the result: 



N^^\t) = GiW{t) (222) 
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FIG. 25: Plot of as a function of G. Fit is to form Eo{G* - G)''^ with values Eq = 1.16, 
G* = 0.79992 andxE = 0.756. 

where 

W{t) = Al f dh r dh{t -h + to)""(ti + to)""(t -t2 + to)'°(t2 + to)~"(ti -t2 + to)"°(223) 
After making the change of variables ti = tx, t2 = ty, e = to/t, then 

W{t) = Alt^-^^'Wie) (224) 

where 

/■I p 1 1 1 1 1 

l^(e) = dx dy- -7 -7 . (225) 

^ ' Jo Jo (1 -x + e)" (x + e)" (1 -y + e)" (y + e)" (x-y + e)" ^ ^ 

It is important to note that the exponent governing the long time dependence is given by 



2 - 5a = 2 - 5(2/3) = -4/3 = -2a and the two terms contributing to N{t) in Eq.(|220|) 
have the same power in time. One can then expect that W{e) is logarithmic in e as e goes 
to zero. A significant amount of work is needed to show that 

W{e) = Wo In (2e)-^ + Wi (226) 
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FIG. 26: Plot of a as a function of G: One-loop case 



where 

Wo = (227) 

" r(2/3) ^ ' 

and the constant Wi could be worked out numerically. The memory kernel is given for long 
times by 

GAl , G,Al 



+ (2e)-^ + ...) (228) 



^4/3 



where 



. (229) 



Exponentiating we have 

{t + toY/^ V 2to 



^W = 7r^f^V(l + ---) (230) 
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We assume that this result will induce a change in the power-law governing the peak ampli- 
tude 

^W = 7^f^r (231) 



and we need to determine u. If N(t) is characterized by the exponent /3 = 4/3 — (5, and 4>{t) 
hj a = 2/3 — u. Following the same steps as at one-loop order we find that the exponents 
satisfy 

(5 - 1 + /3 - 1 = (232) 

gives the result i/ — —S and 

a ^ 2/3 + 5 . (233) 

In evaluating 6 we need the results from the one-loop analysis: — l/(27r-\/3G'), the value 
of Qo with Gi = V2G^/{nQo). One then finds 

Notice that this does no depend on the value of G*. The exponent is increased substantially 
in going from one to two-loop order. More importantly the two-loop theory serves as a 
controlled correction to the one-loop theory. 



F. Numerical analysis: Two-Loop Theory 

We can numerically solve the two-loop amplitude equation rather easily. We expect the 
analytic solution of the last section to hold at the critical point. We first determine the time 

t*{G), when in an unstable run, 0(t) hits a minimum and A^in = (l){t*). We find outstanding 
fits: f* = 1.308/(G'- 0.7032)0-5^1 ^^^^ ^^.^ ^ 0.756(G' - 0.7030)°-^^^ which gives a good first 
estimate for G* — 0.7031 Next we work in the stable phase and compute 

t(G) = / dt (pit) (235) 
Jo 

as a function of G. We find a very good fit to r = 1.482/(0.703235 - Gf-'^^^^ which gives 
an accurate determination of G* — 0.703235. We can then determine (f){t) for G — G*. The 
resulting data can be fit to the form given by 

-Et 

(pa) = Ao-, ^ (236) 
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and we find the outstanding fit with = 0.498, to = 0.206 , a = 0.7377, and E = -0.00019. 
The fit is over the time range <t < 2000. At G — G* we determine 

r(t) = /* dx(f)(x) = Toi^^ (237) 
Jo 

where Tq = 1.040 and = 0.3399. 

The two-loop theory is very similar to the one-loop theory. The analytic work suggests 
a larger shift in the exponent a than is found numerically. One may need to use a self- 
consistent method to obtain more quantitative analytical results. 
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X. KAWASAKI REARRANGEMENT 



A. General Discussion 



We discuss here an approach, due to Kawasaki, which allows one to reinterpret perturba- 
tion theory such that one obtains an ergodic-nonergodic transition at one-loop order. After 
establishing and exploring this result at one-loop order we investigate the stability of this 
solution at two-loop order. 

The kinetic equation for the Laplace transformed correlation function C{z) is given by 
(suppressing the wavenumber dependence in this section) 



where our convention for Laplace transforms is given by 

roc 

CJCU)) = -% / dte'^^C(t) , 

for convolutions we have 

A( /* dsA{t - s)5(s)) = 2/:,(A(t))£,(5(t)) . 



(238) 



(239) 



(240) 



and for time derivatives: 



L,{C{t)) = -i\zC{z)-C\ . 



(241) 



With these results it is easy to see that the inverse Laplace transform of Eq.(238) is given 
by 

Cit) - iK^'^C{t) - f ds K^^\t - s)C{s) = . (242) 
Jo 

Eq.(242) is not of the conventional mode coupling form. KawasakP suggested that the 



kinetic equation, Eq.(238), be rewritten in the form 



z + 



1 + K(')N{z) 



C{z) = C 



(243) 



comparing with Eq.(238) we can solve for N(z) to obtain 

irW(2) 



Niz) 
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If we define 



No(z) 



we can write 



N{z) 



1 _ K(')No{z) 



Eq.(238) can then be written in the form 



(244) 



(245) 



;i + K^'^N{z)){zC{z) -C) + K^'^C{z) = 



(246) 



Taking the inverse Laplace transform gives 



C{t) = -UC(t) -Lo ds N{t - s)C{s) 

Jo 



where 



(247) 



(248) 



sets the time scale. Eq.(247) is of the conventional mode-coupling form. 



We can develop perturbation theory in the dimensionless coupling g as in the development 
above. First we determine K^'^\z) in a power series in g as in previous sections. We insert 



this result in Eq.(245) 



N{z) = [iV(')(z) + iV(')(z) + ...][! + K^'^nI;'\z) + ...] 



N^^\z) + N^^\z) + . . . 



(249) 



with the lowest order approximation given by 



N'^^Xz) = Ni''\z) 



(250) 



and at second order in g we have 



(251) 
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B. One-Loop bare theory 



At one-loop order the mode-coupling kernel is given by 
In the time regime, putting in the wavenumber dependence, 

r(^)(g,t) 



No{q,t) 



where 



1 



{2TTy z + iDr[k^ + {q-ky] 
In bare perturbation theory we have in the structureless approximation 



ir(^'2)(Q,T) = \iVo(g,r) 



where 



(2i) 

In the absence of a cut-off this can be integrated to obtain 



2 \T 



-1 



T{d/2) 



(252) 



(253) 



(254) 



(255) 



(256) 



(257) 



C. One-Loop self- consistent theory 

From the work above we have the self-consistent result at one-loop order: 

if('^'2)(g,r) = g^^iVo(K,T) (258) 

where 

NoiQ, T)=gj ^J{K, T)/(Q - K, T) (259) 
The kinetic equation in this case is given by 



fiQ,T) = -Q'fiQ,T) + Q' dSN,iQ,T-S)fiQ,S) . 

Jo 



(260) 



Consider the difference between the MCT expression, Eq.(260) and the direct solution given 



by Eq.(161). If we use 



/ = -Q'f 



(261) 



valid at lowest order in g, on the right-hand side of Eq.(260) we return to Eq.(161). 
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D. ENE Transition 



The solution for the non-ergodic phase can be separated out as follows. In Laplace 



transform space we have Eq.(243)at one-loop order and in terms of dimensionless variables: 



z + 



1 + iQH~^Nq{Q,z) 



where 



NoiQ,z) 



poo f (j'^pr 

-ij^ dte^^'gJ—fiK,t)fiQ-K,t) 



In the nonergodic phase to leading order for small z, 

HQ) 



and 



NoiQ,z) 



H{Q) 



Inserting these results into Eq.(|262|) and taking the small z limit leads to the result: 

1 



1 + 



and 



H{Q) 



F{Q) 



(27r) 



(262) 



(263) 



(264) 



(265) 



(266) 



(267) 



This set of equations can be solved iteratively. Using comparable numerical methods as 



used to treat the direct approach, we can solve for F(Q) with the results as shown in Fig, 27 



The critical coupling is given by g*^^^ = 9.41. Notice that the wavenumber dependence is 
monotomic. 

E. Two-loop theory 



At two-loop order we need the result 



Nt\z) 



rW(z) 



(268) 
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FIG. 27: Plot of amplitude of nonergodic factor F{Q) for g = 9.41, 9.43, 9.5, 10.0 from the bottom. 



where F*^^^ = F^"* is given by Eq.(157). 

In the nonergodic phase we have the result 



C(Q, u) = F{Q)C{Q)2'k5{uj) + Cr{Q, uo) 



(269) 



where Cr is regular for small u. Inserting this result into Eq.(268) we have as a leading 
contribution for small z: 

r(4) 



NE 



(270) 



where F^^^ is independent of z. This leads to the result 



(4) _ ^4 



q 



(271) 



where is independent of z. Checking order by order we have for the nonergodic phase: 



n2 



(272) 
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where 



N('){z)^K^^^ + ^ . (273) 



Clearly as ^ — > the A^^^-* term dominates the second-order term N^^\ Clearly the ENE 
transition is not a solution at two-loop order. 



65 



XL CONCLUSIONS 

We have introduced a simple dynamic model for a system under going diffusive dynamics 
with a density dependent diffusion coefficient. In the case where the diffusion coefficient 
has constant and linear times in the density, we set up perturbation theory in terms of the 
coefficient of the linear term. For the dynamic structure factor we have worked out the 
associated memory function to fourth order. Analysis of this perturbation theory led us to 
the following conclusions in the simplest case where the static structure factor is a constant 
up to a cutoff: 

(1) . As one increases the dimensionless coupling one finds significant slowing down. The 
observed diffusion coefficient is not modified by higher-order terms in perturbation theory. 
It decreases with increasing density if Di is negative. This gives a mechanism for making 
the coupling g large. 

(2) . For large enough coupling there is a transition where the system goes from stable to 
unstable. 

(3) . Near but below the transition, a slow Fourier component appears that sharpens to 
a 5-function but with an algebraically decaying amplitude. 

(4) . The sharpening of this structural peak corresponds to a new length in the problem 
which grows algebraically with time. 

(5) . Near, but above the transition, the system is metastable with a slow increase with 
time of the amplitude of the peak. Eventually the peak grows exponentially with time and 
the system is rendered unstable. 

(6) . The kinetics of the peak amplitude can be investigated by assuming the peak can 
be approximated by a gaussian with a narrowing width. This leads to a zero-dimensional 
model analogous to the Leutheussar modeP^in MCT. This model can be studied analytically 
near the transition for both one and two loop models. Similarly this model can be studied 
numerically. The emerging picture of power-law decay near the transition is consistent with 
the picture found for the full field theory. 

(7) . We show, for this model, that the ergodic-nonergodic transition, supported at one- 
loop order is not a solution at two-loop order. 

While we have worked out the perturbation theory for a general static structure factor, 
we have explicit results for the simplifying structureless approximation. This corresponds to 
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a coarse-grained model restricted to wavenumbers below the first structure factor peak. The 
resulting kinetic model depends on a single dimensionless parameter the coupling g. At one- 
loop order, as we increased g, we found a critical coupling g* = 93.0 . . . which appears not to 
be a small parameter! However, when we look at two-loop corrections in bare perturbation 
theory for small Q and z, we find a correction, compared to 1, given by C^g and in three 

dimensions C3 — 0.00946 At the critical coupling this gives a correction of 0.45 which 

is acceptable. One explanation for the robustness of perturbation theory is that one could 
introduce the effective coupling g^ff — g/{2T:Y which corresponds to a critical coupling 
^:,, = 0.45.... 

This model is too simple to compare directly with experiment. This is because one 
needs to include the physics at the length scale of the structure factor maximum. One then 
expects an interplay between the mechanism discussed here which controls the generation 
of a metastable structural peak and the peak in the static structure factor. 

The calculation here was carried out in equilibrium. The same techniques can be used 
to treat the associated nonequilibrium quench problem. Also, a similar calculation can be 
carried out for models with density and momentum fields. That case should be interesting 
since the memory function is of the MCT form without rearrangement. Finally this model 
is sufficiently general, diffusion with field dependent diffusion coefficient, that there should 
be additional applications beyond colloids. 

Acknowledgements: This work was supported by the National Science Foundation under 
Contract No. DMR-0099324. 
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APPENDIX A: i2o(^)</'fei</>fc2 • • • 

In developing perturbation theory for time-correlation functions, we need to work out 
the effect of the zeroth-order resolvant operator acting on products of fields. We need to 
evaluate 

W{12...n) = Ro{z)(Pk,K---4>k„ ■ (Al) 
We determine this quantity using the identity 

zW{12 ...n) = <Pk,<Pk, • • • - Ro{z)tDfUk,<Pk, • • • ■ (A2) 
It is not difficult to show that 

n 

DfUk,4>k2 ■■■4>k„=Yl Lo{i)4>kAk2 ---(pk^- Sp (7(12)0^2 ...(j)kj (A3) 



where Lo{l) = Lo(ki) is defined by Eq.(84) and 



7(12) = 2r 'ro(ki, ks) = -2/5-iDoki ■ k2(27r)'^5(ki + ks) 

= (Lo(l) + Lo(2))C(12) , (A4) 
and Sp is an operator which symmetrizes the product it acts on such that 7(ij) appears 



with all possible pairs. Using Eq.(A3) in Eq.(A2) gives 



zWil2 ...n) = <Pk,<Pk, . . . 0fc„ - ^ E Loit)Wil2 ...n) + Sp (^7(12)1^(34 . . . n)) . (A5) 

1=1 

This can be put in the form 

W{12 ...n)= To(12 . . . n) [<i)k,K • • • </>fc„ + (^7(12)W^(34 . . . n))] (A6) 

where 

This allows the W's to be determined recursively. We need W{1) through 1^(1234): 

W{1) = To(l)0(A;i) (A8) 

W{12)=Ul2)[^kA,^-C(kM] + ^^^^ . (A9) 
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W{123) = To(123)[0,,(/)fe,0fe3 - (7(23)0^, - 0(13)0^ - C{12)<f>,,] 

W^(1234) = To{1234:)[cf>k,cf>k,Mk, - {MkJkM] 
+ i?,(1234)] + i^^i^^^^^ 



(AlO) 



(All) 



where 



S,(1234) = S(12, 34) + S(13, 24) + S(14, 23)+ = S(23, 14) + S(24, 13) + S(34, 12|A12) 
and 



B{12, 34) = (To(34) - To(1234))C'(12) 0^30^, - ^(34) 



(A13) 



APPENDIX B: L»JfV(ki)(/'(k2) 



We need to evaluate 



where 



and 



zDj0(ki)0(k2) 



KbT- 



50(x) " 50(x) 



(5 



Ar<^(x,y)^^0(ki)0(k2) 



^;(^)(ki)0(k2) +^;(^)(k2)0(ki) + 5(ki,k2) 



>5(ki,k2) ^-ij d^x J d'^ykeT Ar^{K,y) 



_5 S_ 



(ki)0(k2) 



''Jd'^xJ d'^ykBTDi5(l)(x)5{x - y)2k2 • kie^ 



3i(ki+k2)-x 



= 2i/3-^L'iki • k20(ki + k2) 



■^(ki,k2,k3)0(k3) 



(Bl) 



(B2) 



^(ki, k2, ka) = 2ip-'D,ki ■ k2(27r)'^5(k3 - ki - k2) 



(B3) 



APPENDIX C: INTEGRALS Jo(2/3) AND cjo(2/3) 



Consider the integrals 







and 



Jo(a) 







The second integral can be related to the first via integration by parts: 



Jn(a) 



-(To(2a - 1) 



(2a - 1^ 

We have from Dwight integrals 858.562 and 858.56Jp^ 
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2T{a) 



1 1 

-7 ^ + i—^ 



cos[a^) sm[a^j 



Then 



Jo(a)ao(a) 



T[a)sin{a'K) 



1 



g-if {3a-l 



2a — 1 V{2a — l)r(a) sin[{2a — \)'K)sin[an) 
For the relevant case a = 2/3, we have 

37]-2 1 

Jo(2/3)ao(2/3) = -^(1/3)^(2/3) szn(7r/3).zn(27r/3) " 

Using sm(7r/3) = v^/2, sm(27r/3) = a/3/2, and 

vr 27r 



r(i/3)r(2/3) 



svn(TT 



/3) 



we have finally 



Jo(2/3)ao(2/3) = -2nV3 



70 



1 W. Gotze and L. Sjogren, Rep. Prog. Phys.55, 241 (1992). S. Das, Rev. Mod. Phys. 76, 786 
(2004). 

^ In the standard MSR field-theoretic description perturbation theory is naturally carried out in 
terms of Feynman diagrams with a loop structure. One loop order is second order in Di while 
two-loop order is of fourth order. 

3 K. Kawasaki, Physica A215, 61 (1995). 

^ G. BiroU and J.-P. Bouchaud, Europhys. Lett. 67, 21 (2004) 

5 S. P. Das, G. F. Mazenko, S. Ramaswamy and J. Toner, Phys. Rev. Lett. 54, 118 (1985). 
^ S. P. Das and G. F. Mazenko, Phys. Rev. A34, 2265 (1986). 

R. Schmitz, J. W. Dufty, and P. De, Phys. Rev. Lett. 71, 2069 (1993) 

8 M. Gates and S. Ramaswamy, Phys. Rev. Lett. 96, 135701 (2006) 

9 D. S. Dean, J. Phys. A29 L613 (1996) 

° K. Kawasaki and S. Miyazima, Z. Phys. B103, 423 (1997). 
^ K. Miyasaki and D. R. Reichmann, J. Phys. A38, L343 (2005) 

2 G.De Dominicis and L. Peliti, Phys. Rev. B18, 353 (1978) 

3 A. Andreanov, G. Biroh, and A. Lefevre, J. Stat. Mech., PO2008 (2006). 
^ B. Kim and K. Kawasaki, cond mat/0610588. 

^ G. Mazenko, cond mat/0609591. 

^ G. H. Fredrickson and H. G. Anderson, Phys. Rev. Lett. 53, 1244 (1984). F. Ritort and P. 
SoUich, Adv. Phys. 52, 219 (2003). 

J. P. Garrahan and D. Ghandler, PNAS 100, 9710 (2003); Phys. Rev. Lett. 89 035704 (2002). 
S. Whitelam, L. Berthier, and J. P. Garrahan, Phys. Rev. Lett. 92, 185705 (2005). 
^ S. Whitelam, L. Berthier, and J. P. Garrahan, Phys. Rev. E71, 026128 (2005). 

R. Jack, P. Mayer and P. Sollich, Mappings between Reaction-diffusion and Kinetically Con- 
strained Systems: A + A <^ A and the FA model have upper critical dimensionality dc = 2. 
cond-mat/060152 9l 

This requires a certain amount of coarse graining since S is infinite in the microscopic case. 



71 



22 



An important result in this model is that 



23 G. F. Mazenko, Phys. Rev. A9, 360 (1974). 

24 G. F. Mazenko, S. Ramaswamy and J. Toner, Phys. Rev. Lett. 49, 51 (1982); Phys. Rev.A28, 
1618 (1983). 

25 H. C. Anderson, J. of Phys. Chem. B106, 8326 ( 2002); B107, 10226 (2003) ; B107, 10234 
(2003). 

2^ Due to the one-particle irreduceable property of the memory function, we can replace v ^ 



in Eqs.([l9| and (|20j). 
27 E. Leutheusser, Phys. Rev. A29, 2765 (1984). 

2*^ H. B. Dwight, Tables of Integrals and Other Mathematical Data, 4th ed., MacMilHon, NY 
(1961). 



